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CONTRACTOR REPORT 


INLET TURBULENCE INTENSITY LEVEL AND CROSS-STREAM 
DISTRIBUTION EFFECTS ON THE HEAT TRANSFER 
IN PLANE WALL JETS 

INTRODUCTION 

Numerical simulation has become established, along with experimental and ana- 
lytical investigation, as a useful tool for the analysis and design of engineering equip- 
ment that encounters fluid flow and heat transfer. An important requirement for 
employing the numerical simulation tool is the specification of the appropriate initial 
and boundary conditions to complement the conservation equations that govern the 
momentum and energy transfers of the flow within the domain of interest. 

One requirement that is crucial to obtaining meaningful results is the proper 
specification of the values of the dependent variables at the inlet of the flow domain. 
Earlier contributions by Sturgess et al. [1], Nallasamy and. Chen [2], and Crawford [3] 
have focused on some aspects of the problem. 

Ideally, the values of the dependent variables at the domain inlet should be 
obtained from measurements. A full complement of measured values of the dependent 
variables at the inlet is, however, almost always unavailable. The numerical analyst 
is thus required to make estimates of at least some of these variables at the inlet in 
order to initiate his numerical code. 

The specific flow problem addressed in the present investigation is the effect 
of the turbulence intensity level and its cross- stream distribution at the inlet on the 
heat transfer prediction in a two-dimensional turbulent wall jet. The results obtained 
from numerical simulation are compared with the experimental data of Hishida et al. [4]. 



FORMULATION 


Governing Equations 

In the present work, the essentially parabolic flow situation of a plane jet is 
described using the fully elliptic two-dimensional form of the Navier Stokes equations. 
Thus, the governing conservation equations for the mean turbulent motion and heat 
transfer obtained by applying the Reynolds decomposition and time averaging of the 
instantaneous continuity, momentum, and energy equations are written as follows: 
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The closure of the Reynolds stresses is achieved using the k-e turbulence model. 
Transport equations for k, the turbulence kinetic energy, and e, its dissipation rate, 
are written as follows: 


and 
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where pG^, the turbulence kinetic energy production, is given by 
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Boundary Conditions 

Solid Boundary 

In the present investigation, the near-wall region of the flow domain is not 
resolved numerically. Instead, the standard wall- function approach is used to impose 
wall boundary conditions on the velocity and temperature as well as the turbulence 
kinetic energy and its dissipation rate. The formulation used for the temperature 
wall function is that presented by Hammond [5]. 

Open Boundary 

The need to specify boundary conditions for an open boundary parallel to the 
direction of mean jet motion, as well as for an exit boundary at some downstream 
location, is a direct consequence of the treatment of the flow as elliptic rather than 
parabolic. The open boundary is located sufficiently far away from the slot from 
which the jet issues. This makes it possible to impose zero normal gradient conditions 
for all the variables at the open boundary. An open boundary located ten slot widths 
from the origin was found to be adequate for achieving the desired boundary condi- 
tion for the present investigation. A similar approach was adopted by Amano and 
Brandt [6] and Leschziner and Rodi [7]. 

Inlet Boundary 

One serious limitation associated with most experimental data of fluid flow and 
heat transfer available in the open literature is the lack of a complete set of data for 
all the variables of interest at the inlet plane of the flow domain. The fluid flow and 
heat transfer numerical analyst is thus almost always constrained to specify as best 
he can the values of the flow variables required at the domain inlet in order to 
initiate his code. 
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In the present study, the effect of the intensity level and cross- stream distri- 
bution profile of the turbulence kinetic energy (TKE) at the inlet plane on the heat 
transfer in plane wall jets is demonstrated. Four different values of the TKE level, 
varying through one order of magnitude, were investigated. Also, its cross- stream 
distribution was varied from a uniform profile to various plausible profiles assumed 
below and similar to the one- seventh power law used to specify turbulent velocity 
profiles. 

The Cross-Stream Turbulence Intensity Profile 

It is generally known from experimental data that for wall-bounded flows the 
TKE is maximum at the wall, being the region of the highest shear production and 
minimum at the center where there is no production. A possible profile of the TKE 
distribution in a fully -developed wall-bounded flow is shown schematically in Figure lb. 

The mean value of the TKE is usually specified as a percentage of the kinetic 
energy of the mean flow at inlet. Thus 




( 7 ) 


where 

u = U + u* . (8) 

However, the maximum and minimum values of the TKE for the k(y) profile in Figure 
lb are not known and these can be varied to investigate their effect on the heat 
transfer. 

The one- seventh power law profile for the turbulent velocity in a two-dimensional 
planar channel, as depicted in Figure la, is 
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where 


u = U(8/7) 
max 7 


( 10 ) 


A plausible profile for the TKE, as shown in Figure lb, and complementary to that 
above for the streamwise velocity, is 
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The minimum TKE, k . n , which appears in equation (11) is related to the mean 
TKE, k which is known from equation (7), in conjunction with the total TKE per unit 
depth convected into the flow domain of interest at the inlet plane obtained from 

w/2 

w U k = 2 j dy u k . (13) 

o 


Note that for convenience, the origin for the integration above is taken at the channel 
centerline. Substituting the u(y) and k(y) profiles into the RHS of equation (13) 
and integrating yields 

k min = k C9/C 6+8)] . (14) 

The effects of various k(y) profiles at the inlet plane can now be investigated 
by varying the parameter 0. A value of unity for 0 yields a uniform k(y) profile. 

The corresponding e(y) profile is specified as 


e = C 3/4 k 3/2 L 


(15) 
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where L is an appropriately specified turbulence length scale. 

In fully-developed axisymmetric flows, the corresponding equations for the u(y) 
and k(y) profiles appropriate to the coordinate systems, shown in Figure 2, are as 
follows : 
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where 


e = 


k /k . 
max min 


(19) 


and 


k min = k t 6/ ( g+5)] • (20) 

Exit Boundary 

The exit boundary is located sufficiently far downstream that the zero normal 
gradient condition can be specified for all the dependent variables. 
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NUMERICAL SCHEME 


The governing partial differential equations subject to the appropriate boundary 
conditions, as presented in the previous section, are solved using the finite- volume 
SIMPLE algorithm of Patankar and Spalding [8] and Patankar [9]. A staggered grid 
scheme in which the momentum vector cells are displaced from the scalar cells, as 
outlined in the above references, was adopted. 

The computation was carried out on a 40 W by 10 W flow domain where W is the 
slot width. A 52 x 32 non-uniform grid distribution was used. 

RESULTS AND DISCUSSION 

The experimental data of Hishida et al. [4] have been used to provide the neces- 
sary physical basis for evaluating the numerical predictions. The results obtained from the 
numerical simulations indicate a significant effect of the magnitude and cross- stream 
distribution of the turbulent kinetic energy on the predicted heat transfer in the 
immediate vicinity of the slot. A comparison of the predicted and measured Nusselt 
numbers for an inlet flow Reynolds number of 1700 and TKE values at the inlet of 
between 0.5 and 5 percent of the mean flow kinetic energy is presented in Figure 3. 

The difference in predicted Nusselt number, due to different TKE intensity levels at 
inlet, is restricted to an x/w value of about 15 from the inlet. All the predictions 
collapse onto a single curve for x/w values greater than 15. The predicted Nusselt 
number exhibited a maximum over-prediction of about 30 percent in this region of the 
flow . 

Curves A, B, C, and D in Figure 3 show clearly the significant influence of the 
intensity level of the TKE at the inlet on the predicted heat transfer in the develop- 
ment region of the jet. A similar but less intense effect is observed in Figure 4 
which shows the effect of the cross-stream distribution of the TKE at a single TKE 
intensity level of 1 percent of the inlet mean flow kinetic energy. These predictions 
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demonstrate the importance of the proper specification of the TKE intensity level 
as well as cross-stream distribution in order to obtain correct heat transfer pre- 
dictions especially in the first 15 slot widths from the inlet slot. 

Figures 5 and 6 are similar to Figures 3 and 4 but at almost double the flow 
Reynolds number. The predicted Nusselt numbers in the latter two figures are 
generally higher than those obtained at the lower Reynolds number flows. However, 
the trends observed at the higher Reynolds numbers are the same as were observed 
for the lower Reynolds number flows. 

CONCLUSION 

A numerical study has been undertaken in which the influence of the turbulent 
kinetic energy intensity level and cross- stream distribution at the inlet on the heat 
transfer in a two-dimensional turbulent wall jet was investigated. Both parameters 
were shown to influence significantly the predicted Nusselt number in the developing 
region of the jet, with the zone of influence restricted to about 15 slot widths from 
the point of issue. However, this influence did not extend beyond the x/w = 15 
location and all the predictions were observed to collapse onto a single curve in the 
fully- developed jet region. The predicted Nusselt number in the fully- developed 
region showed a maximum over-prediction of about 30 percent when compared with 
the experimental data. 

A proper specification of the turbulence kinetic energy at the inlet, preferably 
obtained from the experiment, is very important for the correct prediction of the heat 
transfer in the development region of a wall jet. 
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(a) u (y) - PROFILE (b) k (y) - PROFILE 

Figure 2. Streamwise velocity and TKE profiles for fully- developed 

flow in axisymmetric channels. 
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Figure 3. Effect of turbulence intensity at inlet for Re 
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Figure 4. Effect of turbulence intensity cross-stream 
distribution at inlet for Re = 1700. 
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Figure 5. Effect of turbulence intensity at inlet for Re = 3300. 
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Figure 6. Effect of turbulence intensity cross- stream 
distribution at inlet for Re = 3300. 
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